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ABSTRACT 

We find strong evidence for outflows originating in the dense molecular core NGC 6334 I(North): a 
1000 Mq molecular core distinguished by its lack of H ii regions and mid-IR emission. New observations 
were obtained of the SiO (2 — > 1) and (5 -^ 4) lines with the SEST 15-m telescope and the H2 (i^ = 1 ^ 0) 
S(l) line with the ESO 2.2-m telescope. The line profiles of the SiO transitions show broad wings 
extending from —50 to 40 km s~ ^ , and spatial maps of the line wing emission exhibit a bipolar morphology 
with the peaks of the red and blue wing separated by 30". The estimated mass loss rate of the outflow 
is comparable to those for young intermediate to high-mass stars. The near-IR images show eight knots 
of H2 emission. Five of the knots form a linear chain which is displaced from the axis of the SiO outflow; 
these knots may trace shock excited gas along the path of a second, highly coUimated outflow. We 
propose that I(N) is a rare example of a molecular core in an early stage of cluster formation. 

Subject headings: ISM:individual (NGC 6334) - ISM:jets and outflows - ISM:molecules - 
stars:formation 



1. INTRODUCTION 

At a distance of 1.7 Kpc (Neckel 1978), the NGC 6334 
molecular cloud is an 11 pc long filament containing a 
remarkable chain of flve luminous star-forming regions 
(McBreen et al. 1979). The youngest of these star-forming 
regions is thought to be the northernmost of the chain, 
designated NGC 6334 I in the McBreen et al. nomen- 
clature. Evidence for the youth of this region is found 
in the sheer concentration of phenomena thought to ac- 
company the earliest stages of star formation: an array of 
OH, H2O, CH3OH and NH3(3,3) masers (Gaume & Mu- 
tel 1987; Forster & Caswell 1989; Menten & Batrla 1989; 
Kraemer & Jackson 1995), luminous mid-IR sources (Har- 
vey & Gately 1983; Kraemer et al. 1999), an ultracom- 
pact Hii region (DePree et al. 1995), a bipolar outflow 
(Bachiller & Cernicharo 1990; Davis & Eisloeffel 1995; 
Persi et al. 1996), and a cluster of stars detected at 1- 
2 ^m (Tapia, Persi & Roth 1996). 

The young stars of NGC 6334 I are embedded in a mas- 
sive dense molecular core which has been mapped in a 
number of spectral lines and the submillimeter continuum. 
Interestingly, these same maps show that NGC 6334 I is 
part of a twin core system: only 2' north (1 pc at 1.7 
kpc) of NGC 6334 I, a second dense core of roughly equal 
size is apparent. This second core, designated NGC 6334 
I(North) [hereafter: I(N)] by Gezari (1982), has perhaps 
the strongest NH3 lines observed in the sky, yielding a gas 
temperature of 30 K and a volume density in excess of 
10^ cm-3 (Forster et al. 1987; Kuiper et al. 1995). Es- 
timates of the total mass of the cloud range from 1000 to 
3000 Mq (Gezari 1982; Kuiper et al. 1995). However, in 
distinct contrast to source I, NGC 6334 I(N) shows no de- 



tectable H II regions, mid-IR emission, or clusters of 2 /im 
sources (Loughran et al. 1986; EUingsen, Norris & McCul- 
loch 1996; Tapia, Persi & Roth 1996). 

Radio observations have detected an H2O maser, a clus- 
ter of Class I CH3OH masers, and a Class II CH3OH 
maser toward I(N) (Moran & Rodriguez 1980; Kogan & 
Slysh 1998; Walsh et al. 1998). The enigmatic presence of 
masers without associated infrared emission or H 11 regions 
has lead to the conclusions that I(N) is a pre-stellar core 
about to undergo star formation (Kuiper et al. 1995) or 
that I(N) contains a protostar in a pure accretion phase 
(Moran & Rodriguez 1980). In this letter, however, we 
present strong evidence for outflows originating in I(N), 
indicating that the I(N) core already contains at least one 
protostar capable of driving a powerful outflow. 

2. OBSERVATIONS AND DATA REDUCTION 

The near-IR data were obtained over three nights in 
June 1998 with the IRAC2b camera on the ESO 2.2-m 
telescope on La Silla, Chile. The weather was photomet- 
ric and the seeing was 1". The detector was a 256 x 256 
pixel NICM0S3 array with a pixel scale of 0.5". A warm 
Fabry-Perot (FP) interferometer with a spectral resolu- 
tion of X/6X ~ 1000 was used with narrow band fllters 
acting as order sorters. Toward 16 on-source and refer- 
ence positions, we obtained 3 minute integrations at ten 
wavelengths, including three wavelengths bracketing the 
H2 (z^ = 1 — > 0) S(l) line, three bracketing the Br-7 line, 
and four wavelengths selected to measure the continuum 
contribution. Calibration standards were imaged at each 
of the FP settings. To avoid internal reflections from the 
FP, dome flat flelds were obtained for each narrow band 
filter with the FP out of the optical path. The data were 
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reduced with a custom program written in the IDL envi- 
ronment. 

The miUimeter-wave observations were carried out in 
May and June 1998 using the 15-m Swedish-ESO Submil- 
limeter Telescope (SEST). The observations were made 
using SIS receivers with Tgys of 100-200 K and 1440 chan- 
nel acousto-optical spectrometers. The SiO (2 -^ 1) and 
(5 — > 4) lines were mapped simultaneously using the dual 
beam-switching mode with a beam throw of 12' in az- 
imuth. The ^^CO (2 -^ 1) line was subsequently mapped 
in position switching mode with a reference position offset 
by (Aa, A^) = (-20', 20'). In September 1999, ^^CO (2 -^ 
1) spectra were obtained toward three positions; position 
switching was used with an offset of (Aa, A(5) = (30', 0'). 
The FWHP beamwidths are 58" for SiO (2 -> 1) and 23" 
for SiO (5 ^ 4), i^CO (2 -^ 1) and ^^CO (2 -^ 1). The ve- 
locity widths per channel (typical RMS noises per channel) 
are 2.39 (0.03), 0.97 (0.04), 0.95 (0.10) and 0.95 km s-^ 
(0.12 K) for the SiO (2 ^ 1), SiO (5 ^ 4), ^^CO (2 -^ 1), 
and i^CO (2 -^ 1) lines, respectively. Using the CLASS 
package, 1st and 3rd order baselines were subtracted from 
the SiO and CO data, respectively. The temperatures were 
calibrated using the chopper-wheel method and placed on 
a main-beam brightness temperature scale using the beam 
efRciences given in the SEST Handbook ver2.1 (1998). The 
pointing corrections never exceeded 5". 

3. ANALYSIS 

We mapped the SiO (5 — > 4) and (2 — *■ 1) lines in a 
5x5 grid with 20" spacings centered on the I(N) core. 
SiO line emission is detected in both transitions at all 
25 map positions. Broad wings are detected in both the 
(2^1) and (5 — » 4) spectra at velocities ranging from 
Visr — —50 km s~^ to 40 km s~^. In Fig. |], we display 
SiO (2 -^ 1) and (5 -^ 4) spectra toward the positions 
showing the strongest blue and red wing emission. 

The velocity integrated SiO (5 -^ 4) emission of the 
red wing, blue wing, and line core are plotted in Fig ^ 
The wings are integrated over the velocity range in which 
wing emission is apparent in both the SiO (5 -^ 4) and 
(2 -^ 1) lines and is clearly separable from the line core. 
Higher spectral resolution SiO (2 — > 1) and (3 -^ 2) data 
show the wings extending into the -15 to 7.5 km s~^ ve- 
locity interval of the plotted line core emission (Tieftrunk 
& Megeath 1999); however, this low velocity wing compo- 
nent cannot be isolated in the low spectral resolution data 
presented here. 

The SiO line core emission is sharply peaked at a posi- 
tion of ai95o = 17''17"34.3^ ^1950 = -35°42'10", coinci- 
dent to within 10" of the peaks of the dense gas tracers 
HC3N (15 -^ 14), CS (7 -^ 6) and NH3(3,3) (Megeath & 
SoUins 1999; Kraemer & Jackson 1999) and of the peak of 
the 350 - 1100 /im continuum emission (Sandell 1999). 
The red wing peaks 10" to the north of the center of 
I(N) (as defined by the peak core emission); the stronger, 
more spatially extended blue wing peaks 20" to the south- 
east of the center. The spatial offset of the red and blue 
wings is also apparent in the position-velocity diagram 
(Fig. H). The positions of the peaks are not dependent on 
the adopted integration intervals or the choice of contours. 
The detection of broad wings which peak within 20" of the 
center of I(N), and which exhibit a spatial offset between 
the blue and red-shifted emission peaks, is strong evidence 



for a bipolar outflow originating in I(N). 
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Fig. 1. — Spectra of the SiO (2 -^ 1) and SiO (5 -^ 4) emission 
linos toward NGC 6334 I{N). Displayed are the spectra at the peak 
positions of the blue and red wings. 

The column density of SiO in the outflow is determined 
independently for each 5 km s~^ velocity interval from 
Visr = —50 to —15 km s~^ for the blue wing and from 
10 to 40 km s^^ for the red wing. Assuming optically 
thin emission and convolving the SiO (5 — > 4) data to 
the beamsize of the SiO (2 — > 1) data (58"), we calcu- 
late the column density in the J = 2 and J = 5 states 
and estimate a rotational temperature from the ratio of 
these column densities (Coldsmith & Langer 1999). The 
rotational temperatures are ~ 10 iiT for all of the veloc- 
ity intervals. Assuming a constant rotational temperature 
between all of the rotational states, a total SiO column 
density is determined for each velocity interval. The 10 K 
rotational temperature is significantly less than the 30 K 
kinetic temperature of the core (Kuiper et al. 1995), indi- 
cating that the SiO is subthermally excited. A statistical 
equilibrium simulation with a large velocity gradient ap- 
proximation shows that for kinetic temperatures of 30 to 
100 K and volume densities of 10^ cm~^ to 3 x 10^ cm~^, 
the distribution of molecules among the rotational states 
can be well approximated by a contant rotational temper- 
ature with a value of ^ 10 K. With these simulations, we 
estimate that the errors incurred by our assumption of a 
constant rotational temperature are less than 20%. 

The derivation of the outflow properties require an in 
situ measurement of the abundance of SiO relative to H2, 
as the SiO abundance can be enhanced by orders of mag- 
nitude in outflows (Acord, Walmsley & Churchwell 1997; 
Codella, Bachiller & Reipurth 1999). We use the "CO 
(2 -^ 1) spectra to determine the column density of if2- 
Toward the peak position of the blue SiO line wing emis- 
sion, the ^'^CO (2 -^ 1) spectra shows a clear detection of 
the blue line wing. We limit our analysis to the veloc- 
ity range of Visr — —20 to —15 km s~^, over which the 
^■^CO (2 -^ 1) wing has a high signal-to-noise and shows 
no troughs indicative of contaminating emission in the ref- 
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Fig. 2. — The SiO (5 — > 4) line emission. Left panel: the line core emission integrated over Visr = —15 to 7.5 km s~^. The contour levels 
are 3.9 to 17.55 by 1.95 K km s" ^. Middle panel: the dashed contours show the blue-shifted line wing emission integrated over Vj^^ = —50 
to —15 km s~^ and the solid contours show the red-shifted gas integrated over Vi^r = 7.5 to 40 km s~^. The contour levels are 2.5 to 4.5 
by 0.5 K km s~^. The open circles mark the 44 GHz Class I methanol masers (Kogan & Slysh 1998), and the filled square marks the 6 and 
12 GHz Class II methanol maser (Walsh et al. 1998). Right panel: the position-velocity diagram along the the diagonal cut shown in the 
middle panel. The velocity resolution is smoothed to 4.9 km s~^ The contour levels are -0.03, 0.07, 0.14, 0.21, 0.35 to 1.40 by 0.35 K 



erence beam. The ratio of the ^^CO to ^'^CO tempera- 
tures at this position and velocity range indicate that the 
^•^CO is opticaUy thin with an optical depth of 0.1. We 
adopt a kinetic temperature of 30 K (Kuiper et al. 1995) 
and assume LTE, resulting in a ^'^CO column density of 
4.6 X 10^"* cm~^ averaged over a 58" FWHP beamwidth 
and integrated over the —20 to —15 km s""'^ velocity in- 
terval. Assuming a ^^C/^'^C ratio of 59 (Wilson & Rood 
1994) and a [C0]/[H2] abundance of lO"'', the column 
density of H2 is 2.7 x 10^° cm~^. For the identical ve- 
locity interval and beamwidth, the SiO column density is 
2.2 X 10^^ cm^^ and the resulting relative abundance is 
[SiO]/[H2] = 8 X 10"^. This value is similar to that de- 
rived in the outflow of the neighboring NGC 6334 I core 
by Bachiller & Cernicharo (1990). The estimated abun- 
dance decreases with an increase in the adopted kinetic 
temperature: a temperature of 100 K yields [SiO]/[H2] — 
3 X 10~^. A similar analysis of the SiO column density 
in the Gaussian-shaped line cores and a H2 column den- 
sity derived from C^^O measurements (Megeath & Sollins 
1999), results in an SiO relative abundance of 2 x lO"^*' for 
the low-velocity and quiescent gas in I(N), showing that 
the SiO abundance is enhanced in the high velocity gas. 

Table 1 



(z^ = 1 — > 0) S(l) emission (Fig. g). Through a careful 
inspection of the images, we have identified eight knots of 
H2 emission toward I(N). The brightest of the knots is the 
nebula detected in broad band 2.2 /Ltm imaging by Tapia, 
Persi & Roth (1996). The lack of any known UV sources in 
I(N) capable of fluorescing H2 molecules, the distribution 
of the H2 emission in widely separated, compact knots, and 
the absence of continuum emission toward the H2 knots 
suggests that the knots are shock-excited nebulae in one 
or more outflows. The knots are to the northwest of the 
SiO outflow, and they possibly trace the continuation of 
the redshifted lobe of the outflow. However, five of the 
knots fall in a linear chain which is displaced from the 
axis of the bipolar outflow detected in SiO. This geometry 
suggests that these five knots may trace a second, highly 
collimated outflow in I(N). 



Property 


Blue Wing 


Red Wing 


Mass 


4.9 M0 


2.0 Mq 


Momentum 


110 Mq km s-i 


46 Mq km s-^ 


Energy 


3.0 X lO-^s ergs 


1.2 X W^^ ergs 


Mass Loss'' 


1.0 X 10-3 Mq yr-i 


0.4 X 10-3 Mq yr-i 


Luminosity'*'^ 


51 L0 


20 Lq 



"Adopting an age of 5000 years (see text). ''The luminos- 
ity is the mechanical luminosity given by the energy over 
the age. 

Adopting a SiO relative abundance of 8 x 10^^, we es- 
timate a total momentum and energy using the equations 
given in Choi, Evans & Jaffe (1993). To calculate the dy- 
namical timescale, we take the outflow radius to be 15" (or 
0.12 pc), half the distance between the emission peaks of 
the line wings in Fig. 0, and divide by the mass-weighted 
average velocity to obtain a value of 5000 years. The de- 
rived outflow properties are listed in Table 1. These values 
are not corrected by an assumed inclination angle. 

Independent evidence for an outflow in I(N) is shown 
in our continuum subtracted image of the 2.12 fj,m H2 



4. DISCUSSION 

The detection of the bipolar molecular outflow is strong 
evidence for ongoing star formation in NGC 6334 I(N). 
The presence of outflows, and the implication that the core 
is not starless, provides an explanation for the enigmatic 
presence of CH3OH and H2O masers in the core. Indeed, 
the detected masers are coincident with the observed high 
velocity gas (Fig. ||). Since the driving source of the SiO 
outflow has not been directly detected, the source prop- 
erties must be inferred from the properties of the out- 
flow, with the constraint that the source luminosity must 
be less than the total core luminosity of 10'' Lq derived 
from submillimeter observations (Gezari 1982). Shepherd 
& Churchwell (1996) have found a relationship between 
source luminosity vs. mass loss rate for luminosities rang- 
ing from 1 to 10^ Lq. 
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Fig. 3. — Left panel: a /C-band image with contours of HC3N 
(15 — > 14) emission overlaid; the contour levels are 5.19 to 17.3 by 
1.73 K km s-i(Megeath, SoUins & Wilson 1999). The NGC 6334 
I region is centered on the southern HC3N clump; the northern 
clump is the NGC 6334 I(N) core. A dense cluster of 2 fj,m sources 
is apparent toward I. Right panel: the continuum subtracted H2 
(f = 1 — > 0) S(l) emission. Overlayed are the contours of the 
blue (4.0 K km s"!; dashed) and red (2.5 K km s"!; dotted) SiO 
(5 — > 4) line wings and the integrated HC3N (15 — > 14) emission 
(5.19 K km s~^; solid). A complex of H2 emission knots is visible 
toward I. Toward I(N), the positions of the eight H2 emission knots 
are marked. The five knots tracing a possible coUimated outflow are 
marked by the vertical lines. 



loss rate of the NGC 6334 I(N) outflow, 
'^, is characteristic of sources with luminosi- 



The mass 
10-3 Mq yr" 

ties of 10'^ Lq. Due to the unknown inclination of the 
outflow, the actual mass loss rate can be substantially dif- 
ferent; however, it is unlikely that the resulting uncertainty 
in the luminosity is more than a factor of ten. Conse- 
quently, the driving source is probably a deeply embedded 



intermediate to high-mass (proto)star. 

It is likely that I(N) is forming multiple stars. The geom- 
etry of the H2 knots is suggestive of at least one additional 
outflow, driven by a second, embedded source. The detec- 
tion of five faint, highly reddened {H — K > 1.6) near-IR 
point sources and three distinct submillimeter condensa- 
tions within the confines of the I(N) core is further evi- 
dence for multiple sites of star formation (Tapia, Persi & 
Roth 1996; SandeU 1999). Additionally, the mass, den- 
sity and size of the I(N) core are comparable to those of 
cluster-forming cores; Fig. ^ shows that the similar NGC 
6334 I core contains a cluster of more than 90 stars (Tapia, 
Persi & Roth 1996). For these reasons, we propose that 
NGC 6334 I(N) is a rare example of a molecular core in 
the early stages of cluster formation. We predict that I(N) 
will evolve into a region similar to the active I core, with a 
pronounced cluster of 2 /im sources, H 11 regions and lumi- 
nous mid-IR sources. A longstanding question is whether 
star formation in clusters progresses in a sequence, start- 
ing with the low-mass stars and ending with the formation 
of high-mass stars. If the NGC 6334 I(N) core is in the 
process of forming a cluster, then the formation of at least 
one intermediate to high-mass star has preceeded the ap- 
pearance of a dense cluster of low-mass stars. This does 
not necessarily imply that high-mass star formation has 
preceeded the formation of most of the low-mass stars; the 
low-mass stars may be hidden by the high extinction in 
the core. These results do imply that the earliest stages of 
cluster formation may be identified not by bright mid-IR 
emission, H 11 regions or dense clusters of 2 /im sources, 
but by the presence of powerful outflows and masers de- 
tectable at radio wavelengths. Surveys for masers and out- 
flows, particularly in regions without mid-IR sources, is a 
promising means for identifying molecular cores at the on- 
set of cluster formation. 
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